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Abstract It has been proved by S.L.Ziglin [1], for a 
large class of 2-degree-of-freedom (d.o.f) Hamiltonian sys- 
tems, that transverse intersections of the invariant manifolds 
of saddle fixed points imply infinite branching of solutions in 
the complex time plane and the non-existence of a second 
analytic integral of the motion. Here, we review in detail 
our recent results, following a similar approach to show the 
existence of infinitely-sheeted solutions for 2 d.o.f. Hamil- 
tonians which exhibit, upon perturbation, subharmonic bi- 
furcations of resonant tori around an elliptic fixed point [2]. 
Moreover, as shown recently, these Hamiltonian systems are 
non-integrable if their resonant tori form a dense set. These 
results can be extended to the case where the periodic per- 
turbation is not Hamiltonian. 

1 Introduction 

In the last 15 years, there has been active interest 
in the study of integrability (or absence thereof) of 
nonlinear dynamical systems based on the analysis 
of their singularities in the complex time t-plane 
[1-6]. 

Singularity analysis is the study of the behavior of 
the solutions of differential equations around their 
singularities in complex time. While, any anlytic 
system of differential equation is locally integrable, 
the different local patches do not, in general, fit to- 
gether globally. The main idea of the singularity 
analysis is to obtain global information on the in- 
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tegrability of a system through the local analysis of 
the solution in the complex plane. Its origins can 
be found in Kowaleskaya's classical work and on the 
Painleve 's classification of second-order ordinary 
differential equations (see, e.g. [3]). However, for a 
long time, Kowaleskaya's work and Painleve 's the- 
ory were consider interesting, if not old fashioned, 
masterpieces in the theory of special functions and 
little attention was paid to them late 1970's, when 
it was noticed that they were intimately related to 
the theory of solitons. The various Painleve tests 
for ODEs which followed this discovery [3-5] are 
based on the formal existence of Laurent expansions 
for the solutions around the movable singularities of 
the solution in the complex plane. 

According to this approach one seeks to establish 
conditions such that all movable (i.e initial condi- 
tion dependent) singularities of the solutions of the 
equations of motion of the system 

dx 

— =x = f(x,t) x = (xi(t),---,x n (i)) (1.1) 

are isolated poles, i.e that (|1.1|) possesses the so- 
called Painleve property [3-5]. If this is the case, 
then all solutions of (1.1) are meromorphic (single- 
valued and analytic everywhere except at the poles) 
and the system is often integrable, in the sense of 
having global, single-valued integrals of the motion. 

On the other hand, if infinitely multi-valued solu- 
tions are found, one expects that such integrals do 
not exist and the system is called non-integrable. 
The presence of infinitely-sheeted solutions (so- 
called I.S.S property) can be deduced either ana- 
lytically, by showing e.g. that their series expan- 
sions near a singularity contain logarithmic terms, 



or, numerically, by integrating ( |1.1[ ) along contours 
enclosing one or more singularities in t £ C. 

One rigorous approach to the connection between 
non-integrability and the I.S.S property was intro- 
duced, a number of years ago, by Ziglin [1]. He 
showed, under some general assumptions, that 2- 
d.o.f. Hamiltonian systems of the form 



H = Hq(xi,X2,I) + eHi(xi,x 2 ,I, < 
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with Hi 27r-periodic in possessing, for e = 0, a 
closed homoclinic (resp. heteroclinic) orbit, which 
joins one saddle fixed point to itself (resp. two sad- 
dle fixed points), exhibit, for a large class of such 
perturbations and < |e| <C 1, infinitely-sheeted 
solutions. Ziglin's remarkable discovery was that 
he was able to relate directly this I.S.S. property, 
using Mel'nikov's theory to the transversal inter- 
sections of the stable and unstable manifolds of the 
homo(or hetero)clinic orbit. 

Furthermore, Ziglin proved that these perturbed 
Hamiltonians are non-integrable in the sense that 
they cannot have a global, single-valued integral of 
the motion, other than the Hamiltonian itself. 

The purpose of this paper is first to review, 
Ziglin's results on transverse intersections of the in- 
variant manifolds of saddle fixed points and how 
these imply infinite branching of solutions in the 
complex time plane and the non-existence of a sec- 
ond analytic integral of the motion. Then we shall 
prove as reported [2], the I.S.S property for Hamil- 



tonian systems of the form (1.2) around an ellip- 
tic fixed point and thus establish a connection be- 
tween infinite branching of solutions and the non- 
integrability of such systems, as long as they exhibit 
subharmonic bifurcations on a dense set of resonant 
tori [7]. 

Our approach follows closely that of Ziglin, in 
that we use one of Mel'nikov 's theorems to show 
that I.S.S is a direct consequence of subharmonic 
bifurcations (at e ^ 0) of resonant invariant curves 
of the unperturbed (e = 0) system. Moreover, we 
also make the crucial assumption that, for e = 0, 
(1.2) possesses only meromorphic solutions in the 
complex domain. Unlike Ziglin, however, we do not 
assume the presence of a saddle fixed point, whose 
invariant manifolds govern the dynamics. 



In section 2, we present some useful backgraund 
information on the structure and importance of Rie- 
mann surfaces, [14]. 

In section 3, we state Ziglin's theorm on the 
transversal intersection of invariant manifolds in the 
perturbed 2 d.o.f Hamiltonian (|1.2j). For e = 0, 
these invariant manifolds join "smoothly" in a single 
separatrix, or homo (hetero) clinic orbit of the com- 
pletely integrable unperturbed problem, whose slu- 
tions have only poles in complex t. For e / 0, how- 
ever, these manifolds interesct at infinitely many 
points and infinitely branched multi-valued solu- 
tions appear, as predicted by Ziglin, with e.g. log- 
arithmic singularities. 

In section 4, we state our main theorem, as de- 
scribed in [2] , prove that the series expansions near 
a singularity contain logatithmic terms using the 
theory of Abelian integrals and we illustrate our re- 
sults on the example of a driven Duffing oscillator. 

In section 5, we extend and apply our results to 
the case of non-Hamiltonian perturbations. Finally, 
in section 6, we end with some concluding remarks. 



2 Riemann Surfaces 

A Riemann surface X is a connected two- 
dimensional topological manifold with a complex- 
analytic structure on it. The latter implies that 
for each point P E X there is a homeomorphism 
(j) : U — ► V of some neighborhood U 3 P onto 
an open set V £ C, and is defined so that any 
two such homeomorphisms <fi, 4> with U H U 7^ 
are holomorphically compatible, i.e., the mapping 
4>o(j)- 1 : 4>(U DU) — ► cj)(U n U), called a transi- 
tion function, is holomorphic. In what follows, the 
homeomorphism (f> will be referred to as a local pa- 
rameter. Any set <j) i of holomorphically compatible 
local parameters such that the appropriate neigh- 
borhoods U. cover the entire manifold X is called a 
complex atlas of the Riemann surface X. The union 
of the atlases that correspond to the same complex- 
analytic structure on the manifold X, i.e., to the 
same Riemann surface X, is again an atlas. This 
property is violated if the atlases making up a union 
belong to different complex-analytic structures or, 
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which is equivalent, to different, yet tope-logically 
identical Riemann surfaces. 

The simplest examples of a Riemann surface are: 
any open subset of a complex plane C, C itself or 
an extended complex plane E = C U {°°}- The 
two canonical examples for the occurence of multi- 
valuedness are the log function (for each z G C\{0}, 
log(z) are the solution of e w = z) and the function 
z l / q (q G N, q > 2, solution of w q = z). For these 
functions, the points and oo are critical, that is, 
there is no meromorhic continuation around these 
points. However, in any regions 

Dj = | z = re ie ;r > 0, 9 G [a, b], a < b < a + 2vr j 

the function / is single-valued and analytic. To de- 
fine these functions one introduces cuts in the com- 
plex plane, across which the function cannot be con- 
tinued. This approach is not satisfactory since the 
functions are not continuous on dDj. The solution 
to this problem is to extend the domain of defini- 
tion, rather than restricting tha values of the func- 
tion. It is exactly this domain, in which these func- 
tions are single- valued, which forms a sheet of the 
Riemann surface, in the covering-space of C \ {0}. 
For example, each sheet of the surface corresponds 
to a particular branch of log(z). 

Meromorhic functions, i.e., non-constant holo- 
morhic mappings / : X — ► E, constitute a mean- 
ingful object of analysis on Riemann surfaces. The 
local notation f{z) = f o </> -1 (,z) of a meromorhic 
function /, for any local parameters is a mero- 
morhic function of the variable z G 4>(U) in the 
usual sense. The meromorhic function / takes 
evey value c G S the same finite number of times 
(with the multiplicity taken into account). A point 
Pq G / _1 (oo) is said to be pole of the function /. 
In the neighborhood of any Pq G X, a meromorhic 
function / can be represented a convergent Laurent 
series: 

oo 

Cjiz-ZoY, z = cf>(P), z o = 0(P o ) (2.1) 

j=-N 

where <j) is a local parameter, the number N > — oo 
and does not depend on a specific choice of <f>. 



In addition to the notion of a function on a Rie- 
mann surface, we introduce the notion of an Abelian 
differential. An Abelian differential on the Riemann 
surface X is a meromorhic 1-form uj, given on X. 
This implies that we can write uj locally as f{z) dz, 
where f(z) is a meromorphic function of z in its 
domain. For any Abelian differential, the notion 
of a pole and that of a zero are precisely defined, 
along with the notions of multiplicities and that of 
a residue: 

res(o;;Po) = c_i, u(P) = £ Cj (z - z ) j dz (2.2) 

Abelian differentials are usually divided into three 
kinds: holomorphic differentials (first kind), mero- 
morphic differentials with residues equal to zero at 
all singular points (second kind), and meromorhic 
differentials of the general form (third kind). 

Any Abelian differential uj on the Riemann sur- 
face X satisfies the closure condition 

duj = (2.3) 

where da; denotes the total derivative of the 1-form 
uj. Therefore, locally, a primitive function for the 
differential uj always exists and can be defined by 
the formula 

n= f uj (2.4) 

JPo 

for any simply-connected domain on X that in- 
volves (in the case of third-kind differentials) no 
singularities of the differential uj. Formula (|2.4| ) con- 
sidered on the whole surface X, defines, in general, 
a multivalued function called an Abelian integral. 
The division of Abelian differentials into the three 
kinds can be extended to Abelian integrals. Locally, 
Abelian integrals of the first kind are holomorhic 
functions, Abelian integrals of the second kind are 
meromorphic functions, and Abelian integrals of the 
third kind have logarithmic singularities: 

uj = ( 1 h • • •) dz — > = • • • + C-ilnz + • • • 

(2.5) 

We will use expression ( p.5[ ) to prove that the series 
expansions of solutions of ( |1.2[) near a singularity 
contain logarithmic terms, (see Section 4). 
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3 Ziglin's theorem on a 2 d.o.f 
Hamiltonian System 

In this section, we shall state Ziglin's theorem [1] on 
the splitting of separatrices in 2 d.o.f Hamiltonian 
systems. Consider the Hamiltonian 



H = Hq{x 1 ,X2,I) + eHi(xi,X2,I, 



(3.1) 



where (x, y) and (J, (fr) are canonically conjuagate 
pairs of momentum-position and action-angle vari- 
ables respectively. With Ziglin [1], we now make 
the following assumptions about (|3.1| ): 

I Ti is real analytic in some domain of x = (x,y), 
\I — Jo | < /x, |e| < \i and 2ir - periodic in (fr. 



II For e = 0, J = Jo, (3J) has two hyperbolic 
fixed points x + , x_ joined by a doubly asymp- 
totic solution x(t) — * x ± as t — ► ±oo. 



Ill djTC (x(t), Jo) > c > for all t and the solution 



z(t) = (x(t),J ,^(t)), 4>{t) 



dl 



■dt' 
(3.2) 



can be analytically continued to the strip 

n = {0 < Imt < 2vr/A + | 

(where A+ is the positive eigenvalue of the lin- 
earized system about x + ) and has no more than 
a finite number of singular points in n. 

IV 7i(z,e) can be analytically continued for com- 
plex z and 

^<»W), §£(*(«,*,» 

are single valued in n, for all (fro G K, where 
z(t,cfro) denotes the solution z(t)) of (3.2) with 
cp replaced by (fr + (fro- 

Theorem 3.1 (Ziglin [1]) Under the above as- 
sumptions and if cL7ii(z(i, (fro)) has nonzero sum 
of residues \2.3() in II (for at least one (fro ), the sys- 



tem (3.1) possesses multiple-valued solutions I(t) 
Jo + el\ + • • • since 



Ah 



h dt 



for some contour Ten. In fact, since for any given 
(fro, going around T changes the value of I\ by the 
same amount AI\, we conclude that I(t) is infinitely 
branched in the complex t-plane, much like a logi 
function. 

The connection between ([T^) above and the split- 
ting of separatrices comes from Ziglin's proof [1] 
that ([O]) implies that the following integral does 
not vanish identically: 



J(^o 



{JJo,J^i}(z(t),<A + 0o)dt/O (3.4) 



where {•, •} denotes the Poisson bracket, and Jfo, JJi 
are related to the original Hamiltonians by solving 
(3.1) for a (single-valued) J on a constant energy 
surface —J = Hq{x, y) + eH\(x, y, (fr) + ■ ■ ■ . Ziglin 
also proves that ( |3.3[) implies that the Hamiltonian 
system (3.1) does not possess a second analytic in- 
tegral indepedent of H for any sufficiently small 

N/o[i]. 

Theorem 3.1 can be generalized to the non- 
Hamiltonian case of a periodically driven system 

Xi = fi(xi,x 2 ) + egi(xi,x 2 ,t), £ = 1,2 

where gi(x 1 ,x 2 ,t) = gi(xi, x 2 , t + T), gi are an- 
alytic in x\, x 2 and the unperturbed (e = 0) equa- 
tions have single-valued solutions and a smooth 
separatrix joining two fixed saddle points (see [12]). 
The main idea is that if the Mel'nikov integral is 
not identically zero, one can always find an ana- 
lytic function of x%, x 2 which is infinitely-sheeted 
in the complex t-plane. 

Of course, splitting of separatrices does not nec- 
essarily mean the appearance of chaos, since, in 
a non-Hamiltonian system, invariant manifolds of 
saddle points can split, for e ^ 0, without inter- 
secting {J {(fro) 7^ for all (fro)- Splitting does mean, 
however, non-integr ability, in the sense of the ap- 
pearence of I.S.S with the type of infinite multi- 
valuedness one finds in functions with logarithmic 
singularities. 



fyWi(z(t,0o))dt#O (3.3) 
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4 Infinitely Multivalued Solutions 
and Subharmonic Bifurcations 

Let U = D x (I — fj,, I + fi) x Ei 1 be the direct 
product of a domain D C M 2 with coordinates x = 
{x l ,x 2 ) and I be an action variable in the interval 
Iq — H<I<Iq + h with an angular coordinate (ft 
on the circle S . Consider Hamiltonian of the form 

ID 



H{z,e) = Ho(x,I) + eTCi(x,I,(p), z = (x,I,(f>) 

(4.1) 

which is real-analytic in x = (xi,X2),I,4>, 
H (x,I) = F{x) + G(I) for x G D, \I - I \ < n, 
|e| < fj, and is 2ir - periodic in <f>. For e = 0, we have 
the unperturbed system: 



xi = d X2 H (x,I) 
x-2 = -d Xl TC (x,I) 

i = -d (j) n (x,i) = o 

We now make the following assumptions 



(4.2) 



Al. The system ( |4.2[) possesses an elliptic fixed 
point at (0,0) and a family of doubly periodic 
(in t E C) solutions q a (t), which rotate about 
(0,0) in the (xi,x 2 ) plane. Let a S [—1, 1] with 
q- l (t) = (0,0) and q a (t) — ► Q(t) as a — ► 1, 
where (? 1 (t) = Q(t) is a boundary of D. 

A2. Region D of the phase plane of the (uncoupled) 
F(x) system is filled with periodic orbits q a (t), 
whose period T a varies continuously with re- 
spect to the energy 7i(z, 0) = h. Each such 
orbit is a level of F : F(x) = h a provided, for 
energy h > h a , the unperturbed (e = 0) sys- 
tem has a corresponding closed orbit given by 
G^Qi - h a ) = I a . 

A3. The Hamiltonian TC(z, e) can be continued an- 
alytically, but in general, not single-valuedly 
to a domain in complex (z, e) - space. 



A4. The functions 



(4.3) 



with z(t) = (q a (t), I) and z(t, (p) = (q a (t), I, <f)+ 
<j>) are single-valued for every 0€ M. 



For some Iq and sufficiently small \I — Iq\, \x— q°\ 
and |e| the equation of the isoenergy surface, 

H(x,I,<f>,e) = h (4.4a) 

has a single-valued solution for I, 

— I = H(x, <f>, e; h) = Hq(x) + eH\{x, ft; h) H 

(4.4b) 

with 
Hq(x, -H (x)) 



h H {x) = G~ 1 {h-F(x)) 
Hi(x, -H (x),(j)) 



(4.4c) 

n(H (x)) (44d) 

This allows us to go, on the surface ( 4.4a[) (for 
sufficiently small \I — Iq\, \x — q°\ and \e\), from 



system (l4.l|) to the reduced system : 

d X2 H(x,4>,e) x' 2 = -d xl H(x,4>,e) (4.5) 



(where primes denote differentiation with respect to 
<fi). We now consider system ([4.5|) in the extended 



phase space, V = D 1 xE) 1 , which is the direct prod- 
uct of a domain D 1 C D containing the periodic 
orbit q a (t) and S 1 (with angular coordinate <f>) and 
define the Poincare map: 

P e : (x 1 {(j)),x 2 (ct>)) — > (x 1 (0 + 2^),x 2 (0 + 2vr)) 

(4.6) 

on the solutions of (4.5). 

Before we proceed to our main theorem, we need 
to establish an important proposition which is anal- 
ogous to Mel'nikov's second theorem on subhar- 
monic bifurcations, but does not assume the pres- 
ence of a saddle fixed point governing the dynamics 
of the system. 

For fixed h, eqs ( |4.5| ) take the form of a pe- 
riodically perturbed planar system since Ho de- 
pends only on (xi,x 2 ) while Hi has an explicit (f>- 



depedence. The unperturbed Hamiltonian (4.1) has 
the form H.q(x,I) = F(x) + G{I) and the period 
of periodic orbits for this system satisfies the reso- 
nance relationship 



2irm 



nO(/t — h c 



(4.7) 
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for relatively prime integer pairs m, n. In terms of Then, if for some (ft, the sum of the residues at these 
the "new time" (ft, however, the period T a of the H poles 
system is given by 

S=Y, res{d 4> Hi(q a (t j )^)} + (4.9) 



2itm 



n 



and e is small enough: 



Proposition 4.1 (see [2]) Let {Hq, Hi} denote the 
Poisson bracket of Hq,H\ with respect to x\,xi and 
consider a periodic orbit of the unperturbed system 
q a {(ft) with period T a = 27rm/n on the resonant 
closed curve TCq = h a of the unperturbed Poincare 
map Pq, cf. (4-6). Fix h > 0, m, n integer rela- 
tively prime and choose e sufficiently small. Then, 
if the subharmonic Mel 'nikov function 



M m / n {4>) 



2irm 



{H ,H 1 }(q a ((ft-(ft),(ft;h)d(ft 

(4.8) 

has j simple zeros as a function of (ft £ [0, 2irm/n), 
the resonant torus given by (q a ((ft — (ft), (ft) breaks 
into 2k = j/m distinct 2irm-periodic orbits, and 
there are no other 2irm periodic orbits in its neigh- 
borhood. 

We now proceed to our main theorem [2]: 

Theorem 4.1 Consider a real analytic 2 d.o.f 
Hamiltonian system (4-1), which can be reduced to 
the form by virtue of assumptions A1-A4, 



and let q a (t) be a periodic solution of the unper- 
turbed problem, e = 0, with TLq = h a . Assume, 
furthermore, that q a (t) is a meromorphic function 
of t £ C, doubly periodic with real and imaginary 
periods T a , iT' a respectively and can be analytically 
continued in the strip 

L = { t £ C : < Imi <T' a \ 

Finally, suppose that d^Tii{q a {t) , (ft) (where by (ft we 
abbreviate (ft + (ft) has a finite number of poles, in- 
side a closed contour Y C II, where II is the period 
parallelogram 



n 



{ t e C : < Ret < T a , < Imi <T' a \ 



(i) The TCq = h a invariant curve undergoes a sub- 
harmonic bifurcation. 

(ii) The system possesses 0(e) solutions, which are 
infinitely- sheeted in t £ C. 



Remark 4.1 The second result, (ii), follows imme- 
diately from Hamilton's equations and assumption 

flOD : Evaluating the solution I(t) = I +eh(t)-\ 

along a closed contour T C II we have to 0(e): 



AIi 



hdt 



VW(t),<«dt 



r J r 

= -2vriS/0 

Hence, I(t) increases, after every circuit around T, 
by a fixed amount AI(e) such that 



lim 



AI(e) 



7^0 



and the system possesses the I.S.S property. The 
quantity AI\ is an Abelian integral on a Riemann 
surface £ (see section 2) and we note that the 1— 
form u) (third kind), given by cj = d^Tii{q a {t), (ft)dt, 
has a finite number of poles and the Abelian inte- 
gral AIi a logarithmic singularities. Then I(t) is in- 
finitely branched in the complex plane of time much 
like a logt function. 



Remark 4.2 Finally, we can show, using the re- 
sults of [7], that under the above conditions, the re- 
duced system (|4.5| ) is non-integrable. This is done 
by noting that due to (|4.7| ), our system exhibits 
subharmonic bifurcations on a dense (n, m) set of 
invariant tori and hence according to the theorems 
given in [7], it cannot possess a second analytic, 
single-valued integral of the motion. 
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Let us now illustrate Theorem 4.1 on the example 
of a Duffing oscillator, which does not have a fixed 
saddle point and is perturbed by a periodic func- 
tion, which preserves the Hamiltonian formulation 
of the system. More specifically, we consider the 
system of o.d.es 



and 



x\ = x 2 



X2 



-X! - xi 3 + e cos cut (4.10) 



Note that we can always introduce an angle variable 
(f> = uot and a conjugate action variable I, such that 
the Hamiltonian of this system can be written in 
the form (4.4b), i.e 

2 2 4 

H = H + eR\ = — + — + — exi cos <p + Iu 

(4.11) 

It is easy to verify that the unperturbed system 
(e = 0) possesses a family of periodic orbits around 
the elliptic fixed point (0,0), given in terms of the 
Jacobi elliptic functions [10, 11] 

q k (t) =A{ cn(Xt,k),-Xsn(Xt,k)dn(Xt,k)) 

(4.12) 

with A 2 = 2X 2 k 2 and X 2 (2k 2 -l) = -1. Clearly, this 
family of periodic orbits can be analytically contin- 
ued in the strip : 



L 



[t £ C : 



< hat < 2K 



where K, K' are the elliptic integrals of first and 
second kind. Since within any closed contour in 
L, the number of poles is finite, this satisfies the 
corresponding requirement of Theorem 4.1, cf.([1.9|). 

As is well-known [8, 9] subharmonic bifurcations 
occur in this problem when the Mel'nikov integral 



M m/n, 



{iT ,i?i}(g fc (i),t + ^)dt (4.13) 



has simple zeros, with nT a = 2irm. We now com- 
pute the integral ( [4 . 1 3j ) by the method of residues 
and find that for n = 1 



M m l l {t Q 



and for n/1 



■l) m ir 2 m irmK' 

=exp[ j sinwto 



K^/2 - 4k 2 



2K 



M m/n (t ) = V2nuj(2smujto - coswt ) smhK'f3 

(4.14) 



S 



E 



res(cn(At,-, k) sin 4>) ^ 



where II is the period parallelogram of the cn func- 
tion. Thus, according to Theorem 4.1, AI\ = 
—2niS ^ and the infinitely-sheeted property of 
the solution I(t) = Io+eIi(t)+- ■ •, to O(e), has been 
established. We recall that the existence of such 
infinitely-sheeted solutions for a Duffing oscillator 
similar to ( 4. lCj ) has been explicitly demonstrated 
elsewhere [12, 13], in the form of so-called psi series 
expansions, involving logarithmic singularities. 

5 The case of non— Hamiltonian 
perturbations 

Consider a Hamiltonian system of the form ([4,l|), 
( [1.2| ) perturbed by dissipative terms as follows: 

±i = d X2 F{x) + ed X2 Hi + £71/1 

X2 = -d Xl F(x) - ed Xx Hi + £72/2 

j> = fi(J) + ediHx + e5 l9l 

I = -e{d ( i > H 1 -S 2 g 2 } (5.1) 

where F, G and Hi are as defined in section 4 and 
fi, gi are analytic functions of (xi,X2,F) and 2ir- 
periodic in (j). Furthermore we assume that the fi, 
gi are such that the above perturbation is of non- 
Hamiltonian character. 

In this case, the function TL = Tto(xi, x 2 , 1) + 
eHi(xi,X2, 1, 4>), with Ho = F(xi,x 2 ) + G(I), is no 
longer conserved, as it satisfies 

T~t = £(lifid Xl 'H + 72f2d X2 H + 52g2diH + 6igid < j ) H) 

(5.2) 

or to order e, 

H = eh( Xl ,X2,I) + 0(.e 2 ) (5.3) 

with h(xi,x 2 ,I) = r yifid Xl F+j 2 f2d X2 F+52g2^(I)- 
We now apply the implicit function theorem and 
solve the equation H(x, I,(p) = H for I 

/ = H (xi,X2;H) + eHi(xi,X2,(f);H) H 

whence, after some calculations, the reduced system 
has the form: 

9H _ £ \]h_ _ 7i i 1 _ 9H 5igr 
8x2 ^8x2 ^ dx2 ^ - 



+ 0(e 2 
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x 2 



dxi 



+ e 



dHi _ 72/2 _ dH 5igx 



Applying this approach to non-Hamiltonian per- 



H = ett 



0292-^ 7l/l 



dH 



+ 0(e 2 ) turbations (|5.l| ) of the example of section 4, we find 
that the formula ([T^) yields: 



72/ 2 J+0(^) (5.4) Mrn/n = _ Xj 



dx\ 8x2 

Let d = (71, 72, #i, £2) denote the dissipation co- 
efficients in our system. Also, let P tt d denote the 
Poincare map associated with the system ( |5.4| ). As 
before, we consider a function TL = H(xi, x 2 , 1, (ft, e) 
which can be continued analytically, but not, in gen- 
eral, single-valuedly to a domain in complex space, 
while the function 



4nK 



sn 2 (u,k)dn 2 (u,k)du (5.10) 



R(q a (t),I,^ + 



{-£-5 2 g 2 }(q a (t),I,cft + 



(5.5) 

is single-valued for every (ft. 

We now extend the main theorem of section 4 to 



(5.1), assuming that the function (|5.5| ) has a finite 
number of poles inside a closed contour T C II, such 
that 

£ res(^ - 5 29 2) + (5.6) 



S 



Thus, we conclude that the non-Hamiltonian sys- 
tem (5.1) possesses infinitely sheeted solutions to 
0(e), of the form I(t) = Iq + eli(i) + • • •, since 



h dt 



- S 2 g 2 (q a {t), } dt = -2iriS ± (5.7) 

Finally, it is easy to verify, using the ideas of The- 
orem 4.1, that condition ( |5.6| ) also implies that the 
Mel'nikov function of system (|5.4|) 



M™ /n {4>) = M m ' n {(ft) + M™! n ((ft) (5.8) 



with 



M™ /n {4>) = 7, 



n 



2-Km 



dxo 



72 j 2 



dHp 
dx\ 



7i/i d<P 



(5.9) 

cf. ( |5.8| ) is not identically zero, for general per- 
turbations fi,gi- Hence, if ( |5.8|) turns out to have 
simple zeros, the perturbed system will exhibit sub- 
harmonic bifurcations as expected. 



where we have taken f\ = 0, f 2 = x 2 , g\ = g 2 = 
cf. ( 4.12| ). Performing the integrations in ( 5.1C| ) we 
finally find 



jQrm/n 

72 ~ l2 Zk 2 ^/i=W 



[E(k)(l-2k 2 )-k' 2 K(k)} 
(5.11) 



When added to ( (L14 ), the constant M™J n for 



n 



1, obtained above for suitable choices of 72 



72 < R m (to) 



m/l 



12 



(— l) m 72\/27rw exp[- 



mirK' 
2K 



(5.12) 



can prevent the total M™^((ft, 72) of ( |5.8| ) from hav- 
ing simple zeros, in </> = u>to, thus eliminate the oc- 
currence of subharmonic bifurcations in this exam- 



ple and the unperturbed periodic orbit ( 4.12 ) satis- 
fying the resonance relation 

2irm 



AK(k)Vl - 2k 2 

6 Conclusions 



In this paper, we have shown that 2 d.o.f Hamil- 
tonians which exhibit, upon perturbation, subhar- 
monic bifurcations of their resonant invariant curves 
around an elliptic fixed point possess solutions 
which are infinitely-sheeted in the complex domain. 
Such systems are known to be non-integrable when 
these bifurcating tori form a dense set. 

Our analysis follows Ziglin's approach to simi- 
lar systems, with the important difference that, in 
this case, the unperturbed Hamiltonian possesses 
a homoclinic orbit. In our systems, we do not re- 
quire the existence of such an orbit and make cru- 
cial use of Mel'nikov's subharmonic function. We 
have applied our results to a conservative Duffing 
oscillator, showing that this classical example ex- 
hibits infinitely-sheeted solutions near an elliptic 
fixed point. 
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Finally, we have shown that our results can be 
extended to the case of non-Hamiltonian perturba- 
tions, and have illustrated the analysis by applying 
to a class of non-Hamiltonian perturbation of Duff- 
ing's oscillator. 
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